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Abstract 

The off-equilibrium response function X^Pi^w) and autocorrelation function 

C(t,t w ) of an Ising chain with spin-exchange dynamics are studied numerically 

and compared with the same quantities in the case of spin- flip dynamics. It 

is found that, even though these quantities are separately different in the two 

cases, the parametric plot of x(i, t w ) versus C(t,t w ) is the same. While this 

result could be expected in higher dimensionality, where x(C) is related to 

the equilibrium state, it is far from trivial in the one dimensional case where 

this relation does not hold. The origin of the universality of x(C) is traced 

back to the optimization of domains position with respect to the perturbing 

external field. This mechanism is investigated resorting to models with a 

single domain moving in a random environment. 
64.75. +g, 05.40.-a, 05.50.+q, 05.70.Ln 
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I. INTRODUCTION 



The generalization of the fluctuation-dissipation theorem (FDT) to slowly relaxing sys- 
tems, such as glasses, is an issue of foremost importance for understanding non-equilibrium 
processes. Ordinary FDT relates the autocorrelation C(t,t w ) and the integrated response 
x(t, t w ) functions, which both depend in equilibrium on the time difference t — t w , where t w 
is the time elapsed after the sample preparation. A general feature of slow kinetics, instead, 
is the aging property, namely the dependence of the time-scale of relaxation on the time t w . 
This feature generally shows up both in C(t,t w ) and • In the context of mean-field 

models for spin glasses it was shown [|TJ that, for large t w , x(t, t w ) depends on the two times 
through the autocorrelation function alone 



This property holds quite generally in a wide class of aging systems where deviations from 
the ordinary equilibrium FDT, namely T\{C) = C(t,t) — C(t,t w ), result in a nontrivial 
fluctuation dissipation ratio X{C) = —Tdx{C)/dC. Recently, a theorem has been proven 
linking X(C) to static properties 



where P{q) is the equilibrium probability distribution of the overlaps. This opens the way 
to a classification of aging systems according to the structure of their equilibrium states |§ 
and to the recognition that the properties of the response to a perturbation are universal in 
systems sharing the same overlap distribution. 

The ordering process of ferromagnetic systems provides a simplified framework for the 
study of the off-equilibrium FDT because the main features of slow relaxation are fully 
exhibited but the structure of their equilibrium state is simple and exactly known. While (|]) 
is generally obeyed [0-0], the validity of (||) in this case depends on dimensionality. For d > 1, 
(0) holds asymptotically, implying the same fluctuation dissipation ratio, and hence of x{C), 



x{t,t w ) = x[C(t,t w )]. 



(1) 



dX(C) 
dC 



C=q 



P(q) 



(2) 
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for all these systems ||. This applies, in particular, to different dynamical realizations of 
the same Hamiltonian model such as ferromagnets with non-conserved (NCOP) or conserved 
(COP) order parameter [|]]. However the picture is totally different in the case of the Ising 
chain. With NCOP, (0) is not obeyed and, instead, a non trivial X{C) is found that 
cannot be connected to static properties. 

These results suggest that the nature of X(C) is essentially dynamical in this case. A 
natural question, then, is about which properties of the kinetics are reflected by x(C). In 
order to address this point, we study in this paper the response of the one-dimensional 
Ising model with Kawasaki spin-exchange dynamics quenched to a low temperature, in the 
scaling regime preceding equilibration. In this case the order parameter is conserved and the 
microscopic mechanism whereby coarsening of domains is produced differs completely from 
the dynamics with single spin flip. With NCOP interfaces are independent brownian walkers 
whose density is progressively reduced due to annihilation events. With COP, instead, 
the motion of interfaces is mediated by evaporation, diffusion and recondensation of single 
monomers. Despite this completely different character of the dynamics, we show that x(C) 
is the same for COP and NCOP. Due to the violation of the hypotheses of theorem (0) 
this universal character cannot be traced back to statics but it is more likely to have a 
common dynamical origin. Since the basic coarsening mechanisms with COP or NCOP are 
profoundly different, other kinetic properties, of a more general and fundamental nature, 
determine x{C) ■ The analysis carried out in this article shows that the total response 
can be viewed as due to the elementary contributions given by single domains. In the 
kinetic process domains coarsen and translate and the complex interplay between these two 
mechanisms produces the response. Starting from this idea we introduce simple models 
where a single domain is allowed to diffuse in a random environment. In this framework, the 
elementary response generated by the domain can be studied and from its knowledge the 
behavior of the original Ising chain is inferred by means of scaling arguments. The details of 
the rules for the motion do not change the overall behavior of the response. This approach 
provides a clear physical interpretation of how the response is produced in the Ising model. 
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These considerations are of a general nature and apply, in principle, to any dimension. 
On the other hand, the occurrence in d — 1 of an off-equilibrium response which never 
vanishes, as opposed to the cases with d > 1, is due to the special character of domain walls 
motion in one dimension. Actually, in d > 1 the evolution of a domain is the result of two 
competing drives. The first is the tendency to lower surface tension by making interfaces 
straight. The second is the drift towards regions where the random field is favorable. For 
long times the first mechanism always prevails, and the response generated by the drift of 
domains is negligible ||. In d — 1, instead, domain walls are point-like and surface tension 
does not play any role; moreover the drift mechanism is so efficient as to generate a non- 
vanishing response even in the limit of large times, when the interface density decreases to 
zero. 

This paper is organized in six Sections. Section [lj is devoted to a description of the 



COP dynamics of the unperturbed Ising chain. In Section |TTT| the effects of a perturbation 
are discussed and the response function is introduced, showing the analogy with the NCOP 
dynamics and the universality of x{C). Models for a single diffusing domain are discussed in 
Sections |TV] and [V], where scaling arguments are presented to illustrate the common origin 
of the response within the two types of dynamics. In Section [VI] we discuss the relevance of 
our results for different systems and draw some conclusions. 

II. UNPERTURBED DYNAMICS 

We consider the one-dimensional Ising model with ferromagnetic nearest-neighbor cou- 
pling constant, whose Hamiltonian is 

N 

7-Co({si}) = ~J^2 s * s m, ( 3 ) 
i=i 

where = ±1. The system is quenched from an uncorrelated high temperature equilibrium 
state to the final temperature T. Evolution takes place through Kawasaki spin-exchange 
dynamics, i.e. swaps between antiparallel nearest-neighbor spins. In this way magnetization 
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is a conserved quantity. The model describes lattice gases or binary alloys. The probability 
of exchanging Si,s i+ i is assumed to be 



Pt = mm 



e T , 1 



(4) 



where T is measured in units of the Boltzmann constant and AE = 2J(sj_iSj + Sj+iSj+2) is 
the energy change. 

The basic features of the dynamics following an instantaneous quench are discussed in 0. 
Depending on the energy change elementary moves can be distinguished into three classes 
with AE = 4J, 0, —4 J. The first kind of process is evaporation, namely the separation of a 
spin from the boundary of a domain. Processes with AE = are the diffusion of a single 
spin (monomer) in the bulk of a domain of the other phase. When a diffusing monomer 
reaches an interface a condensation event occurs: The spin joins a domain. This process 
implies an energy change AE = — 4J. 

Evaporation is an activated process occurring over a characteristic time r ev = exp(4J/T). 
For low temperatures r ev is large and one observes a long interval t <C r ev during which 
evaporation practically does not happen. In this regime a reduction of the kink density 
p(t) can be obtained only by the diffusion and condensation of the monomers present in 
the initial state. In order to do this single spins move a distance of the order of the initial 
coherence length £ in a typical time £ 2 . This leads to a decay |7||| of the density of diffusing 
single spins over a characteristic time t s ~ £ 2 . In the case we are concerned with, a quench 
from very high temperature, t s is of order unity and one observes a fast decay to a plateau 
on short time-scales (Fig. p. In the regime r s < f < r ev no diffusing spins are left and 
pit) remains constant. At times of order r ev evaporation events begin to occur and the 
dynamics restarts (Fig. [I]). Evaporated spins diffuse and they may recondense on a kink 
different from the one where they were emitted. This is the well known mechanism leading 
to the decay pit) oc {t/T ev )~ l l z 0: in this regime, dynamic scaling is obeyed. This behavior 
lasts until t ~ t^ op such that p(t^ op ) = p eq ~ exp (— 2J/T), the equilibrium kink density. 
At this time domains reach a size such that a second monomer is emitted when the first 



one is still diffusing. When they meet they form a stable dimer and this process exactly 
balances the domain annihilation due to the evaporation-condensation mechanism, so that 
p keeps its equilibrium value. Regarding the data presented in Fig. [I], for T = 0.7 the 
exponent 1/3 is not clearly observed because the system equilibrates too soon. For T = 0.48 
the effective exponent gradually decreases toward 1/3. At the longest times simulated the 
effective exponent is ~ 0.3. 

The autocorrelation function is defined as 

1 N 

C{t,t w ) = —J2(si{t)si(t w )), (5) 

i=l 

where (■ ■ •) indicates thermal averaging. This quantity is shown in Fig. [2|. C(t,t w ) strongly 
depends on the range of times (t w ,t) considered. In the case t w -C t^ op considered here the 
behavior of the autocorrelation function is different for t -C T^ q OP or t 3> t^ op . For t -C T^ q OP 
C(t,t w ) decays as [p(t) / p(t w )] x . The exponent A depends on t w as follows: When t w = 0, 
the upper bound A < d, originally proposed by Fisher and Huse |J provides the correct 
value A = 1 for the conserved d = 1 Ising model, as shown analytically and numerically 
TU| . Instead, for t w chosen well inside the scaling regime, Yeung, Rao and Desai JTTJ 



in 



found a lower bound A > 3/2 for d = 1. Since A = 1 for t w = this constraint implies the 
dependence of A on t w . To our knowledge, there are no results for the actual value of this 
exponent when t w ^ 0. From the data presented in Fig. |] one observes a power-law decay 
consistent with C(t,t w ) ~ (t w /t) 1 / 2 for the cases t w = 10 4 , 10 5 . Recalling that p{t) ~ i~ - 3 in 
the range of time considered, one obtains A ~ 1.66 > 3/2. In this way we show that the lower 



bound determined in |TTJ is correct and that the value of A with t w chosen inside the scaling 
regime is different from the case with t w = 0. Actually, the value A ~ 1.66 may indicate 
that the value A = 3/2 could be asymptotically correct. In order to check this point lower 
temperatures and larger waiting times should be considered. By plotting C(t,t w ) against 
p{t)/p{t w ) we have checked that also the curve with t w = 10 3 gives the same exponent 
A > 3/2, whereas the smaller exponent observed in Fig. |2] is simply due to p(t) decaying 
with an effective exponent considerably smaller than 1/3 in the range of times plotted in the 
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Figure. For t w = 10 6 the curve starts decaying with the same exponent A = 3/2 but then 
the decrease becomes faster, indicating that the system is close to reaching equilibrium. A 
second important observation, regarding Figure 0, is the convergence of C(t, t w ) towards the 



scaling behavior C(t,t w ) = C(t w /t), as expected quite generally for slow relaxation |L2] and 
in particular for coarsening systems [|T^]. However, differently from NCOP, the convergence 
in this case is very slow and very large t w must be considered in order to exhibit a good 
data-collapse. 

III. RESPONSE TO A PERTURBATION 

Let us consider the Ising model quenched to temperature T in zero field. At time t w a 
random field 

hi = hei (6) 
is applied, so that the hamiltonian is changed into 

N 

n = n -J2 h iSi- (7) 

i=i 

The field takes randomly only two values, £j = ±1, with expectations 

el = (8) 



6i€ 



kr (9) 



The probability of exchanging two spins Sj, s i+ i is ((D with AE = 2 J(sj_iSj + s i+ iSj +2 ) + 
2h(si€i + Si + \€i + i). We consider h/T sufficiently small in order to be in the linear response 
regime. We are interested in the scaling regime, i. e. times such that p(t) ^> p eq . Moreover 
we want the qualitative features of the dynamics, presented above, to be unaffected by the 
external field. This imposes an additional constraint p(t) ^> ^ 1 (h), where = AJ 2 /h 2 
is the Imry-Ma length |T^]. For the values of h and T considered, p eq so that 
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p(t) is unchanged by the presence of the random field from the instant of the quench up to 
equilibration (Fig. [I]). 

We consider the integrated response function 



N 



X(t,t w ) = lim ^77-J2 e i( S i)hi ( 10 ) 

where (■ ■ -)h denotes average in presence of the external field. Before discussing the response 
of the model with COP, let us briefly recall the behavior with NCOP. In this case equilibrium 
is reached because spins are flipped spontaneously in the bulk of ordered domains; the 
characteristic time for this process is T^ q COP = exp(4 J/T). The asymptotic value of the 
response function is the equilibrium susceptibility 

Xe Q = l/T, (11) 

namely lim^oo Tx(t, t w ) = 1. In equilibrium the bulk of domains is responsible for the 
response because spins anti-aligned with the random field are more likely to be reversed by 
thermal excitations. However, for times t < spins in the bulk do not flip; the response 

observed in this regime is then of non-equilibrium nature. Since the bulk is frozen, x(t, t w ) 
must be provided by interfaces. Specifically the motion of kinks is such as to optimize the 
position of domains with respect to the random field, building up a finite response ||. The 
range of times over which the non-equilibrium pattern is observed can be expanded by letting 
T — > (keeping h/T small, for linear response theory to hold), or J — > oo. We refer to the 
latter limit for simplicity, that has the advantage of being easily implemented numerically 
by forbidding the flip of spins in the bulk. With J = oo, t^ cop = oo and the system never 
equilibrates. The exact solution of the model H with J = oo yields the aging form 



7T V t 

that converges, in the large t limit, to 



V% ft 

Tx(t,t w ) = arctanW 1 (12) 



T Xoo = ^= (13) 



S 



for large t. For finite J, as already anticipated, the response of the model is the same as 
with J = oo for times t <C T^ q COP while, for larger times, the equilibrium susceptibility is 
recovered. 

For what concerns the generalization of the FDT, notice that ( |3~2| ) obeys the scaling 
form x(t,t w ) = x(t w /t). Hence eliminating t w /t with C(t,t w ) = C{t w /t) with J = oo one 
finds 



x(C) = — arctan 

7T 



v2cot f |c 



(14) 



This curve is plotted in Fig. [| In the limit C(t,t w ) — > 0, namely t — ► oo, the value Xoo 
is recovered, as previously discussed. For finite J the behavior of C(t,t w ) and xi^^w) is 
unchanged with respect to the case J = oo up to t^ cop '. Therefore, in the fluctuation- 
dissipation plot in Fig. ^| the same curve as with J = oo is followed from C(t w ,t w ) = 1 
down to C eq = C{r^ q COP , t w ). For C(t, t w ) < C eq , namely for t > t^ oop , the system goes to 
equilibrium, x{C) departs from the master curve with J = oo and approaches Xeq- For fixed 
t w , C eq grows with temperature, because r^ q COP decreases as T is increased. Alternatively, 
for a given temperature, C eq grows by increasing t w . In conclusion, the master curve is 
followed in a wider range by decreasing t w or T, because in this way C eq is reduced. This 
explains the behavior of x(C) in Fig. [3[ 

For Kawasaki dynamics the kinetic process is more complex than in the non-conserved 
case. The different coarsening mechanism does not simply change the growth law exponent, 
but even the two-time quantities considered here are radically modified with respect to 
NCOP. In particular, as discussed in the previous Section, the exponent A = 3/2 differs 
from the value [0 A = 1 for Glauber dynamics. Given these differences the fact that the 
fluctuation-dissipation plot turns out to be the same for NCOP and COP, as will be shown 
below, is unexpected and far from being trivial. 

The fluctuation-dissipation plot with COP is shown in Fig. Curves with different t w 
collapse on the same master curve for C(t, t w ) > C eq . The collapse of curves with different t w 
proves the validity of (H) in the scaling regime even with COP. Moreover, this master curve 
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is the same as with NCOP. As for NCOP, the collapse occurs only when both times t and 
t w belong to the scaling regime. With NCOP scaling is obeyed starting from a microscopic 
time to, which is temperature independent [[HJ; with COP this regime is entered after r ev , 
which depends on T. Therefore, for low temperatures the time scales over which curves 
with different t w coincide are completely different in the two revealed from the 

waiting times reported in Figs. ^| and 0. Comparing the two Figures, one also concludes 
that changing t w or T only produces a shift of C eq , the point where x(C) deviates from the 
master curve, in complete analogy for both dynamics. Then, on the basis of what is known 
with NCOP, we expect in the zero temperature limit the master-curve to be followed down 
to C(t,t w ) = 0. Let us stress that the correlation and the response functions are different 
for the two dynamics. Only when the response x{t-> tw) is expressed as a function of the 
correlation C(t,t w ) one obtains the same fluctuation-dissipation relation xiP)- 

The universality of the non-equilibrium response with respect to the type of dynamics 
naturally raises the question of a possible common fundamental origin. As proposed in 
Section Q, this must be of a dynamical character, since the connection (|2p between statics 
and x{C) cannot be invoked in d = 1. The main dynamical feature of phase-ordering is the 
presence of a coarsening structure with many competing domains. This suggests to look for 
the underlying universality in the response of a single domain to the perturbation. In order 
to test this idea, in the following Sections we investigate simplified models for the motion 
of a single domain. Underlying this approach is the assumption that x(t>tw) °f the whole 
system can be seen as the sum of the response of single domains considered independent. 
Correlations between domains are only responsible for the growth of their typical size. In 
this way we are able to identify the diffusive wandering of domains as the origin of the 
non-equilibrium response, both for NCOP and COP dynamics. 
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IV. A SINGLE DOMAIN MODEL WITH RIGID DIFFUSION 



Let us consider an isolated domain T>i of up spins covering the segment [j, j + I] of an 
infinite one- dimensional lattice. A quenched random variable hi is defined on the sites of 
the lattice via (HHH). At each time step T>i is allowed to move rigidly one lattice unit on the 
right or on the left with probability given by (U), with AE = hj — hj + i +1 or AE = hj + i — hj_i 
respectively. The model can be regarded as an Ising chain with an initial condition containing 
I up spins in the interval [j, j + /] in a sea of down spins. This Ising model is governed by a 
dynamical rule which conserves the magnetization and allows only rigid translations of the 
up domain. AE is then exactly the energy gain computed through the Ising Hamiltonian (^) 
with i running over the sites occupied by the domain. 

For h = the landscape is flat and the position x of the center of the domain performs 
a random walk. Provided the linear response regime (h/T — * 0) is considered, also for finite 
h the root mean square displacement Ax(t,t w ) in a time interval [t w ,t] obeys Ax(t,t w ) ~ 
y/t — t w , as shown in the inset of Fig. [|. The response function of the domain is defined as 



X Vl (t 1 t w )= lim -(]>>>,. (15) 

5,->0 II ieDi 

The notation (■ ■ -)h indicates averaging, for a single realization of the random field, over the 
trajectories of T>\. Notice that, differently from (|10D, i runs only over sites of T>\. 

Let us consider the behavior of 'X Dl (t : t w ) for short times. Initially, at t = t w , the response 
is zero. By moving one lattice spacing the sum in ( |ToD can change by a value 0, +2, — 2. 
In the first two cases the energy is unchanged or decreased; in the last one AE > and 
the probability that such a move is accepted is exp(— 2h/T). Then, averaging over the 



trajectories and the random field, (|T5|) gives Tx Vl (t w + l,t w ) = lim.h^ [T /h](J2 i&Vl ei)h = 
limh_^ T[l — exp(— 2h/T)]/(2h). Taking the limit h/T — > the linear response function 
follows Tx Vl (t w + l,t w ) = 1. 

For long times X^K^?^) approaches, assuming equilibration, the static susceptibility 



x ^ = IZh dh (16) 
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where 

■ ! . h(l-2k) 

Z(h,T) = Y,Pke-— (17) 

is the partition function. Here k is the number of sites inside the domain where = 1 and 
Pk = is the probability of having a particular value of k. From fllTD one easily finds 



Z{h,T) 



cosh , 



(18) 



yielding T% v ' = \im h/T ^ (lT/h)t a nh(h/T) [cosh (/i/T)]'" 1 = I. 



The behavior of y^ l {t, t w ), obtained numerically, is plotted in Fig. |, showing the validity 
of the scaling form 

Tx V '(t,t w ) = lg(y) (19) 

where y(t,t w ) = Ax(t,t w )/l. The scaling function behaves as 

{y for y 1 
(20) 
1 for y > 1 

in agreement with the analytical results for short and long times. 

So far we have studied a model where the size I of the domain is conserved by the 
dynamical rule. However, in order to apply this result to the description of the Ising chain, 
where domains coarsen, we consider now a slightly modified version where the size of the 
single domain varies stochastically while growing on average. In such a situation, we indicate 
with l(t) the size of the domain at time t, with L(t) the average of l(t), and define y{t,t w ) 
via 

y(t,t w ) = Ax(t,t V} )/L(t). (21) 

We have performed numerical simulations where the size l(t) of T>\ was increased with a 
stochastic rule such that L(t) = L(0) + a^/i, with L(0) = 10 and a = 1/8 (these values 
are chosen for numerical convenience). Fig. |5| shows that even in this case the scaling 
form fll9y20|) holds. 
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A. Connection with the Ising model 



From the knowledge of the response of a single domain, we can recover the behavior of the 
Ising chain, where many domains compete, by assuming that the latter can be adequately 
described by a collection of quasi-independent domains of average size L(t). This means 
that all effects produced by correlations between domains, apart from the increase of L(t), 
are supposed not to be relevant. This assumption will be further discussed in Section |V|. 

The overall response of the Ising model is then given by 

^g^p(M)f(t,g (22) 
i 

where V(l, t) is the fraction of spins belonging to domains of size I at time t, which obeys [18 



the scaling form V(l,t) = L(t) _1 /[//L(£)]. The analysis can be carried out more easily with 
the approximation f(x) ~ 5{x — 1), i. e. assuming that all domains have exactly the same 
size L(t). Then 

X (t, t w ) ~ L{t)- l x VL {t, t w ) ~ -g{y) (23) 

where y(t,t w ) is defined by ( j21~|) and Ax(t, t w ) is the average distance traveled by domains 
of size L(t). On the basis of the large time behavior of y(t,t w ), three situations can be 
distinguished, namely lim^oo y(t, t w ) = oo, lim^oo y(t, t w ) = const. > or lim^oo y(t, t w ) = 
0, giving rise to different values for the response. 

The first case occurs in the equilibrium state of the Ising model because domains diffuse 
while their size stays constant due to the formation of new kinks. On the basis of (P^D, for 
the total response of the Ising chain one obtains lim^oo x(t, t w ) = Xeq = 1/T. This is indeed 
the value ( |Tl~D found for the original Ising chain. 

The second case, lim^oo y(t, t w ) = const. > 0, occurs when Ax(t, t w ) and L(t) grow with 
the same exponent. This happens in the scaling regime, both with spin-flip or Kawasaki 
dynamics. With NCOP this is true because interfaces are brownian walkers 

Ax(t,t w ) ~ (t-t tt )3, (24) 
13 



while it is well known Jl3| that 

L(t) ~ t*. (25) 

With spin-exchange dynamics the behavior of Ax(t, t w ) can be obtained with the follow- 
ing argument: The displacement of domains is mediated by the evaporation and condensa- 
tion of monomers. When a monomer travels a distance L(t), leaving one domain and joining 
the nearest, one of the boundaries of each of the two domains is displaced by one lattice 
unit. The number of monomers per unit time that leave a domain and reach the neighbor 
instead of recondensing on the original one is p(t) 0. Hence, on average, it takes a time 
~ L(t) to move a domain by a unitary distance. Since domains move randomly the law of 
brownian motion 

[Ax{t,t w )f = 2D(t-t w ), (26) 
is obeyed, with a diffusion coefficient 

D~ p (t). (27) 
Using the appropriate growth law for COP 

L(t) = p^y 1 ~ P (28) 

one gets 

Ax(t,t w ) ~ [{t-t w )/t*y. (29) 

Hence for long times Ax(t,t w ) ~ tP, with = 1/2 for NCOP and (3 = 1/3 for COP. 
Comparing ( p4| , p5| ) and ( p8| , p9"D one concludes, as already anticipated, that Ax(t,t w ) ~ L(t) 
regardless of the dynamical rule. Hence, recalling ( pl|) one obtains with both dynamics the 
universal form 

V = Voc(l-^)* (30) 
14 



where is a constant that may depend on the dynamical rule. Inserting this expression 
into (p3| ) one recognizes that the global response obeys the scaling behavior x(t, t w ) = x(t w /t) 
which, as discussed in Section is correct for the Ising chain. Moreover one also obtains 
lim^ooy = i/oo < oo and, therefore, an asymptotic non-equilibrium value Xoo = g(yoo)/T ^ 
Xeq is generated. This limiting value is different from the static susceptibility Xeq because 
2/oc is finite, as opposite to the equilibrium case discussed above where lim^oot/ = oo. This 
is in agreement with the behavior of the original Ising model (|13|), with both dynamics. 

Our approach does not allow the evaluation of y^. However it reproduces the main 
features of the non-equilibrium response and offers an insight into what goes on after a 
perturbation has been switched on in the one- dimensional Ising model. In particular, the 
model points out clearly which is the mechanism whereby the response is produced. Actually, 
the domain responds to the perturbation by moving so as to optimize its position with respect 
to the random field. Such a sharp statement is made possible in this context by the fact 
that the dynamical rules do not allow any other possibility and indicates that the same 
mechanism is at work also in the Ising chain. The discussion presented insofar shows also 
that another physical ingredient plays a fundamental role in d = 1: The convergence to a 
finite value of y(t,t w ), namely lim^oo y(t, t w ) = y^ ^ 0. This property holds because the 
displacement of domains in d = 1 is proportional to their average size. With these ingredients 
the one-domain model predicts a finite limiting value lim^oo %(t, t w ) = Xoo of the response 
of the Ising chain. We emphasize that the property lim^oo y(t, t w ) = y^ ^ is far from 
being trivial. Although we restrict the analysis in this paper to one dimension, we believe 
this to hold only in the d = 1 case. As discussed in Section |, in one dimension domains 
diffuse in order to lower the magnetic energy in absence of the additional force produced by 
surface tension, because kinks are point-like objects. In d > 1, instead, the displacement of 
an interface is not only ruled by the random field but is also governed by curvature. This 
additional mechanism lowers surface tension and competes with the tendency to lower the 
magnetic energy. The weakening of the drift of domains limits their motion so that Ax(t, t w ) 
grows more slowly than L{t). As a result lim^oo y(t, t w ) = and the response produced in 
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this way, from ([19]) and (|20|), vanishes. Therefore we expect the third possible behavior of 
y(t,t w ) introduced above to be realized for d > 1. 

V. A SINGLE DOMAIN MODEL WITH SPIN EXCHANGE DYNAMICS 

In the previous Section we have discussed the diffusion of a rod of average size L(t). 



With this dynamics we have obtained the formula for the response. The connection 
with the Ising model was then made possible by scaling arguments where Ax(t, t w ) and L(t) 
were assumed a priori to behave as in the NCOP or COP Ising model. The dynamics of the 
rod however is quite different from the actual behavior of the Ising model, where generally 
domains do not move as a whole. Furthermore, while diffusion is the mechanism whereby 
interfaces move in the Ising model with NCOP, as stated by fl24|) , with COP the law (|29|) 
is obeyed, showing a non brownian character. Despite these shortcomings of the model, we 
have obtained a good description of the Ising chain and this suggests that the actual details 
of the dynamics are irrelevant. 

In this Section we introduce an improved single-domain model which takes into ac- 
count the different kinetics for NCOP and COP. In particular, for COP, the evaporation- 
condensation mechanism which rules the evolution is taken into account. 

We study an Ising chain of size iV with periodic boundary conditions and two interfaces 
initially located in x\ — 1, x% — N/2 + 1. We consider both spin-flip and spin-exchange 
dynamics, and require the initial structure with only two domains of opposite sign to be 
preserved at all times. This condition is guaranteed if the temperature is low enough. 
However simulations at very low T would be numerically too demanding, in particular for 
COP. For this reason we implement the aforementioned condition in a different way: With 
NCOP we let J = oo by forbidding flips in the bulk. For COP we use a modified dynamics 
where new dimers are not allowed to form. Specifically, when two monomers meet, one of 
the two particles is removed and attached randomly to another kink. With these rules a 
configuration with only two domains persists. 
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Let us focus on one of the two domains, indicated with T>, with the center located in 
(xi + x%)j1. With COP the system evolves via exchanges of monomers between the two 
boundaries of T>. Since at low T at most one monomer is present the size I = X2 — x\ of T> 
is practically conserved. With NCOP, on the other hand, interfaces diffuse independently 
and the size changes. We consider the range of times t <C (N/2) 2 so that annihilation events 
do not happen. Although the size of domains changes, the average value L = (x2 — x\)h is 
constant, due to symmetry. Then L = const, with both types of dynamics. Our goal is to 
compute 




Notice that this is exactly the response (|10D defined for the Ising model. The subscript T> 
simply reminds that we are in a situation with only two domains of fixed size. 

Before discussing the behavior of the model let us comment on the relationship between 
this approach and the one presented in the previous Section. An obvious difference is the 
presence of two domains instead of one. However, with the choice L = N/2 the domains 
are equivalent and this merely doubles the response. On the other hand, in the regime t <C 
(N/2) 2 considered now, y(t,t w ) = Ax(t,t w )/ L <C 1. Therefore, if this model is equivalent to 
the previous one we expect to recover the results with y(t,t w ) <C 1 of the previous Section. 

With NCOP the behavior of the model can be deduced from the knowledge of the re- 
sponse Xsing(t, t w ) of the case with a single kink located in x(t) and fixed boundary conditions 
discussed in ||. For this system it was shown exactly that 

Y n 2 

Xsing(t,t w ) = Jim -T77-J2( S i)h e i = l^ Sx (^ *w) ~ (t - t w )* (32) 
A^oiV/i- =1 Nl 

where 5x(t,t w ) = (| x(t) — x(t w ) \) is the average distance traveled by the kink in the 

time interval Result fl32|) allows one to deduce the behavior of the present model. 

Denoting by xiin 9 (^^) an d xiin 9 (^^) the responses associated to the two interfaces the 

total response is simply given by 

X V (t, t w ) ~ x2U*, tw) + xfiLgit, tw) ~ Sxi(t, t w ) + 5x 2 (t, t w ) (33) 
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because the two interfaces are independent. Indicating with Ax(t,t w ) = [8xi(t,t w ) + 
8x 2 (t,t w )]/2 oc (t - t w f/ 2 the average distance traveled by V in the interval [t w ,t] one 
finds 



X v {t,t w ) ~Ax(t,t w ) ~ (t-t w )*. (34) 

Going back to the Ising model, assuming again that domains are non-interacting, the 
response is obtained by multiplying x V (tytw) times the number of domains present p(t) ~ 

r 1 / 2 , 

X {t,t w ) ~p(t) X V (t,t w ), (35) 

yielding 

X(t,t w ) = Xoo (l - y) ' ( 36 ) 



where a = 1/2 and Xoo is the asymptotic value. Then, from (|T^,|^,^) one recovers the 
behavior of the previous model in the small y(t,t w ) limit, as expected. 

With NCOP this result has been obtained by letting the single interface move as in the 
original Ising model, namely with the same update rules for the spins. With COP this issue 
is more subtle. In the Ising model the diffusivity of domains depends on time via (p?|). In 
the present case, with a fixed size of the domains, D is constant. To keep this into account, 
we consider a spin-exchange dynamics generated by the modified probability 
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p T (s l} s i+ i) = min \n s(t)e T , lj , (37) 

where n(t) is a counter of evaporation events. With this rule the diffusivity is proportional 
to n- 1 / 2 ^) ~ p(t), as in the Ising model. The response of the model, obtained by numerical 
simulations, is plotted in Fig. and for long times 

X V (t,t w ) ~ (t-t w )K (38) 

From ( |3"8"| ) the response of the Ising model is obtained through (P5|), leading to the same 
form (^) as for NCOP, but with a = 1/3. This shows that the present model, evolving with 



two different dynamical rules with NCOP or COP, gives rise to different responses. This 
agrees with the behavior of the Ising model, where both an d C(t,t w ) are different 

but x(C) is the same. 



The approach in terms of single domains is based on the assumption (|35|) of their quasi- 
independence. We have shown that this hypothesis allows a description of the Ising kinetics 
in terms of scaling arguments and provides a good agreement with the original model. We 
are now in a position to substantiate further the validity of this assumption by checking the 



accuracy of (]35f) . 

From the knowledge of x{t,t w ) for the Ising model we extract the effective response due 
to a single domain defined by 

X(t, t w ) = p(t)xeff(t, t w ). (39) 

The accuracy of the independent domain approximation can be determined by comparing 
X T> (t,tw) with Xeff(t,t w ). With NCOP this issue has been considered in || showing a very 
good agreement. In particular x V {t,t w ) and Xeff(t,t w ) both increase as t 1 / 2 for large t. For 
COP the behavior of Xeff(t,t w ) is shown in Fig. [7|. For large t, in particular, y^{t,t w ) an d 
Xeff(t,t w ) grow with the same exponent 1/3. 

The results of this Section point out the robustness of the mechanism generating the 
response in d — 1 which only relies on the coarsening domain structure of the system. Pro- 
vided this character is maintained the global behavior of the susceptibility and, in particular, 
the convergence to a finite non-equilibrium value, is the same. 



VI. DISCUSSION 

In this paper we have studied the off-equilibrium response of the Id Ising model. We have 
shown that the fluctuation dissipation plot is the same with NCOP or COP. In d — 1, where 
the connection (0) with statics cannot be invoked, this universal character has a dynamical 
origin, as shown by the analysis of simplified models presented in Sections |TV] and |V|. 
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An important issue is the relevance of the picture provided by the Ising chain for arbitrary 
dimension. In d = 1 with NCOP the Ising model reaches equilibration on the characteristic 
time t^L = p~ q 2 = exp(4 J/T). As mentioned in Section [TT1] this is the time necessary 
for flipping spins in the bulk of domains; the same mechanism is also responsible for the 
equilibrium response Xeqit — t w ). The off-equilibrium response, instead, develops in the 
regime t < r^ q COP . Hence in d = 1 two kinds of response exist which are observed on 
different time-scales separated by t^ cop . This feature gives rise to the pattern presented 
in Fig. §. 

The case d > 1 presents some differences. For quenches below T c global equilibration 
is never reached in an infinite system. Despite this fact the response can still be split into 
an equilibrium and an aging part: The bulk of domains, which behaves as a pure phase 
and attains local equilibrium, produces Xeq(t — t w ) which obeys (0) and (0). Domain walls, 
instead, are responsible for the non-equilibrium part which obeys the scaling form 

Txit,t w )=r w a f(^] (40) 



til 



with 



^ , for d < 3 

(41) 

| , for d > 3 

and logarithmic corrections in d = 3. The dependence of a on dimensionality results from 
the competition between the drift of interfaces produced by the perturbation and the force 
caused by their curvature. When a = 1/2, as for d > 3, the response is simply proportional 
to p(t) implying that a single interface produces a response which does not depend on 
time. This is what happens if interfacial spins simply polarize according to the random field 
on a microscopic time-scale. On the other hand, from the knowledge of the behavior of 
the one-dimensional case, we know that the wandering of interfaces gives rise to a single- 
interface response growing as (t — t w ) l l 2 . Therefore a natural interpretation is the following: 
Curvature, which is absent in d = 1 becomes progressively more important as d increases, 
due to the coordination number. The attempt to lower surface tension weakens the drift 
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of domain walls and inhibits the response mechanism associated with it. This progressively 
increases the value of a with respect to d — 1 as dimensionality is increased. Then, for 
d > 3 the motion of domain walls is fully governed by curvature whereas, for d < 3 the drift 
mechanism partly compensate the decrease of p(t), resulting in a smaller exponent a. Only 
at the lower critical dimension d = 1, however, this mechanism is so efficient as to balance 
exactly the loss of interfaces yielding a = and an asymptotic finite limit Xoo- 

For d > 1 the presence of an equilibrium and an off-equilibrium response, and the mech- 
anisms whereby they are produced, strongly resemble the situation in d = 1. However, while 
for d = 1 they are observed on different time-scales, for d > 1 they are both developed during 
the phase-ordering process. Then, for d > 1, since the equilibrium part alone obeys (0), in 
order for the total response to fit into the scheme (|2]), the off-equilibrium contribution must 
vanish in the large t w limit. ( [41]) shows that this happens for d > 1 but the decay of the 
off-equilibrium response for d < 3 is slower than usually expected on the basis of the idea 
that the random field simply polarizes the interfacial spins if d < 3. 

For d = 1 in the phase-separation regime Xeq{t — t w ) is absent and only the off-equilibrium 
response is developed which in this case does not vanish for t w — > oo and causes the break- 
down of the connection (^j with the statics. 

In conclusion, with NCOP an overall discussion of the response of the Ising model to 
stochastic perturbations can be given in terms of two mechanisms whose interplay is regu- 
lated by dimensionality. The equilibrium response, which prevails in d > 1, only relies on the 
structure of the equilibrium state through (0) and, therefore, is independent from dynamics. 
In this paper it was shown that in d = 1 also the off-equilibrium response is independent 
on the kinetic rules, although this property has a different origin. The possible universality 
of the out of equilibrium response in higher dimensionality and the generality of the scaling 
form (f4"0|) are interesting issues which deserve to be investigated in the future. 

As a final comment, we discuss the possible relevance of our studies for systems with 
a vector order parameter with N components. In this case spins basically rotate rather 
than flip and this is totally different with respect to scalar systems. The absence of bulk 
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and interfaces prevents a straightforward extension of the concepts developed in this paper 
and the mechanisms by which the response is built up in the vectorial case is complex and 
still not well understood. However, the exact computation of the response function in the 
solvable large- iV model |HJ has recently shown a pattern which resembles the behavior of 
scalar systems. Actually in the large- N model the response function can be explicitly split 
into an equilibrium and an off-equilibrium part. It can be shown, moreover, that the former 
satisfies (||D while the latter obeys ( flUl) with an exponent a that vanishes at the lower critical 
dimension of the model which, due to the vectorial character, is d = 2. This close analogy 
with the scalar case shows that, although the microscopic dynamics is different, the same 
competition between two mechanisms for the development of the response exists probably 
for every value of N and that the same scaling relation ( fiO| ) holds. This suggests the idea 
that the off-equilibrium response may be independent from the dynamics also in vectorial 
systems, as we have shown for the Ising model at the lower critical dimension. 
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FIGURES 
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FIG. 1. The kink density p is plotted against time for quenches with T = 0.48 and h = (black 
solid line) and T = 0.7, both with the perturbation (h = 10~ 2 , black bold solid line) and without 
(h = 0, white solid line collapsing on the curve with h = 10~ 2 ). Data refer to numerical simulations 
of a system of N = 10 5 spins, averaged over 10 realizations. The dashed line represents the i -1 / 3 
law. 
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FIG. 4. The integrated response function is plotted against the autocorrelation function with 
COP. Data are obtained from numerical simulations of a system of N = 10 5 spins with h = 1CP 2 , 
averaged over 4000 realizations. The dashed line is the case with NCOP and J = oo, as in Fig. |^. 
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FIG. 5. The scaled single domain response Tl 1 x Dl (t,t w ) is plotted against Ax(t,t w )/l for 
T = 0.7, h = 10~ 2 and t w = 0. Circles, squares and diamonds correspond to domains of size 
I = 10, 20, 40, respectively. Averages are taken over 10 6 trajectories. The dashed line is the 
analytic behavior for small Ax. Triangles are the response of a domain whose average size L(t) 
grows according to L(t) = 10 + yi/8. In the inset the behavior of Ax(t, t w ) is plotted for a domain 
of size I — 40 with t w — 0. The dashed line represents the law (i — t w ^j ' . 




FIG. 6. The response of a pair of domains diffusing via Kawasaki dynamics with probabil- 
ity (|37|). Simulations are presented for a system of iV = 10 5 spins with T = 0.7 and h = 10~ 2 
averaged over 2 • 10 5 realizations. The dashed line is the i 1 / 3 behavior. 
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FIG. 7. The effective response of the Ising model with Kawasaki dynamics for a quench to 
T = 0.48 with t w = 3 ■ 10 4 and h = 1(T 2 . Data are averaged over 4000 realizations. The dashed 
line is the t 1//3 behavior. 
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